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FLUCTUATING PROPERTIES OF TURBULENT BOUNDARY LAYERS 
FOR MACH NUMBERS UP TO 9 


By William D. Harvey, Dennis M. Bushnell, 
and Ivan E. Beckwith 
Langley Research Center 

SUMMARY 

Estimates of the intensity of the longitudinal velocity fluctuations across a hyper- 
sonic turbulent boundary layer have been made on the basis of experimental measure- 
ments of fluctuating density and pitot pressure. These data were obtained at the Cornell 
Aeronautical Laboratory in the boundary layer on the wall of a shock-tunnel nozzle at 
Mach numbers from 8.2 to 8.9. Both the fluctuating and mean density data were obtained 
from the intensity of light emitted by the interaction of an electron beam with the local 
airflow. 

A review of previous hot-wire data and assessment of the new data indicate that 
reasonable estimates of both the magnitude and trends of the intensity of fluctuating 
density, velocity, and temperature can be made by application of a simple mixing-length 
concept to mean flow data. Comparison of the new results with the previous hot-wire 
data indicates that the root-mean-square values of these fluctuating quantities may be 
essentially independent of Mach number up to about Mach 9, provided that suitable nor- 
malizing parameters are used. 


INTRODUCTION 

Most of the available experimental data on compressible turbulent boundary layers 
is limited primarily to the mean properties of the flow. In order to obtain a better under- 
standing of such boundary layers and ultimately advance theoretical prediction methods, 
it is necessary to study the statistical properties of the flow. 

Measurements of mass flow and total temperature fluctuations in supersonic turbu- 
lent shear layers have been made with hot wires (see refs. 1 to 3). Root-mean-square 
levels of density, static temperature, and velocity fluctuations were obtained from these 
hot-wire data by applying a small-perturbation analysis and neglecting static pressure 
fluctuations. Direct measurements of density fluctuations and mean density have been 
made in hypersonic turbulent wakes by the use of electron-beam (ref. 4) and interfero- 
metric (ref. 5) techniques. Recent measurements of the mean and fluctuating density and 



pitot pressure across a turbulent boundary layer on a shock-tunnel wall at a nominal 

Mach number of 8.5 were made at Cornell Aeronautical Laboratory and have been \ 

ij; 

reported in reference 6. The density and pitot -pressure fluctuation data used herein j 

were obtained from the original analog records for the Cornell Aeronautical Laboratory j 

research. The purpose of the present paper is to estimate the intensity of the velocity j 

fluctuations from these data of reference 6. Also an evaluation of a mixing-length j 

approach to compute the magnitude and trend of the density and velocity fluctuations j 

from data of reference 6 and other compressible data is made. Finally, comparisons 
are made with previous fluctuation measurements in zero-pressure-gradient boundary 
layers to determine whether any Mach number or Reynolds number effects can be 
established. 

SYMBOLS 

A second-order density-velocity correlation coefficient defined in equation (8) 

a = (y - 1)M 2 



Cp specific heat at constant pressure 

G function of Mach number defined in equation (5) 

H total enthalpy 

h static enthalpy 

l mixing length 

M Mach number 

n exponent in equation (B2) 



Pt 


2 


r 


p 

Pt 


pressure 
pitot pressure 


q = (u 2 + 2uu’ + q’ 2 ) 1//2 


q' = (u' 2 + v'2 + w' 2 ) 1 / 2 


R gas constant 


^u 


second-order temperature-velocity correlation coefficient, 


R o,R o ,R 9 ,R 9 third-order correlation coefficients defined 

uq^ pu“ } pq^ p“u 

Rg Reynolds number based on boundary-layer thickness 


T'u* 

CTu.2 

in equations (C2) 



u 

u,v,w 


distance along streamline 
temperature 

local stagnation temperature 

stagnation point downstream of normal shock 


longitudinal, normal, and transverse velocity components, respectively 



x 


y 

OL 

y 


longitudinal coordinate 
normal distance from tunnel wall 

angle between local fluctuating velocity vector and pitot tube axis 
specific heat ratio 


3 



6 


nominal boundary-layer thickness 


fji viscosity 

p density 

r shear stress 

Subscripts: 

aw adiabatic wall 

e mean values at edge of boundary layer 

o settling chamber conditions 

s edge of sublayer as denoted by peak in \/u'2 or in static values 

w wall conditions 

a pitot -tube angle of attack 

1,2 conditions ahead of and behind normal shock, respectively 

A bar over a symbol indicates a time mean value. A prime denotes a fluctuating 
quantity. 


APPARATUS AND EXPERIMENTS 
Tunnel and Test Conditions 

Measurements of density and pitot pressure were made in the turbulent boundary 
layer on the nozzle wall of the Cornell Aeronautical Laboratory 48-inch hypersonic shock 
tunnel (the actual diameter of the nozzle at the data station was about 23.5 inches or 
0.597 meter). Data were obtained at a nominal Mach number of 8.5, and measurements 
of skin friction, heat transfer, and pressure on the nozzle wall were made simultaneously 
(ref. 6). All boundary-layer-profile measurements were made at a single station 
135 inches (3.429 meters) downstream of the nozzle throat. A summary of the test con- 
ditions for the experimental program is provided in reference 6 where the free-stream 
conditions in the test section were calculated by assuming equilibrium air expansion to 


4 


r 


the Mach number corresponding to the measured test-section pitot pressure. For con- 
venience, these test conditions from reference 6 for the four runs used herein are listed 
in the table which follows (in the notation of the present paper). The assumption of 
equilibrium expansion can be justified on the basis of data reported in reference 7. 


T 

o 


Po 

Pw 

Pt 

M e 

T 

e 

Free -stream 
Reynolds 
number 
per foot 
(per 30.5 cm) 

°r 

°K 

psia 

MN/m2 

psia 

N/m^ 

psia 

N/m2 

°R 

OK 

5805 

3225 

845 

5.8 

*0.035 

*241 

3.30 

22 752 

8.25 

502 

279 

1.52 X 10 5 

3795 

2108 

240 

1.6 

*.0115 

*79 

1.10 

7 584 

8.62 

272 

151 

1.10 

2920 

1622 

290 

2.0 

*0125 

*86 

1.29 

8 894 

8.89 

189 

105 

2.12 

3590 

1994 

385 

2.6 

*016 

*110 

1.61 

11 100 

8.82 

245 

136 

1.84 


Instrumentation 

The mean and fluctuating density measurements were obtained with a 40-kilovolt 
electron-beam density probe which was projected through an orifice in the nozzle wall. 

A graphite cup was used to collect the beam current of approximately 1 milliampere. 
Photomultipliers were used to detect the fluorescent light emitted from collisions 
between the gas molecules and the beam electrons. A discussion of the design, opera- 
tion, and calibration of the electron beam is included in the appendix of reference 6. 

The mean and fluctuating pitot pressures were obtained by using piezoelectric 
crystals mounted in 0.125-inch-diameter (0.003175-meter) pitot tubes which were cov- 
ered with protective caps. There was a small volume of gas between the crystals and 
the caps which had a small orifice to transmit the stagnation-pressure pulse to the crys- 
tals. A more complete discussion of the pitot probes used may be found in reference 6. 
Several possible types of interference between probes or between the viewing optics port 
and the wall instrumentation were evaluated, and the results of the evaluation may also 
be found in reference 6. The recording system, which was used to obtain the pitot data, 
sampled a given signal every 50 microseconds, stored it on a drum, and afterwards dis- 
played the signal on a strip-chart recorder. 

CALCULATION OF FLUCTUATING QUANTITIES 


Mixing-Length Equations 

A mixing-length formulation has been used to model temperature and velocity fluc- 
tuations for incompressible flow (refs. 8 and 9, for example) and velocity fluctuations 
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for compressible flow (ref. 10). An analogous mixing-length relation for the root-mean- 
square fluctuating density may be written as 


l 

1 "o 
| a> 

A (p/Pe) 

p 5 p 

d(y/ 5) 


Similarly, the root-mean-square fluctuating temperature may be expressed as 


]/t' 2 1 

d(T/T e ) 

T e 5 

d(y/6) 


and the root-mean-square fluctuating velocity, as 


_ l 

d(u/u e ) 

u e 6 

d(y/ 5) 


( 1 ) 


( 2 ) 


(3) 


where for the present purposes Prandtl's mixing-length concept is used in the inner por- 
tion of the boundary layer - that is, 


l = 0.4y 



Recent calculations (refs. 9 and 10) indicate that the distribution for l in the outer part 
of the boundary layer can be approximated as 


l = 0.086 



Slopes of the profiles of mean density, temperature, or velocity are used herein with 
equations (1), (2), or (3) to compute fluctuating values. 


Expressions for Velocity Fluctuations 

The extension to compressible flow of recent calculation methods for incompress- 
ible flows that utilize the turbulent-kinetic-energy equation (refs. 11, 12, and 13) is to a 
large extent dependent upon a knowledge of the fluctuating quantities for compressible 
flow. Both the magnitude and distribution of the root-mean-square fluctuating velocities 
in a turbulent boundary layer are important parameters required in such methods. 

A general relation (see appendix A) between the instantaneous velocity vector, den- 
sity, and pitot pressure for values of a up to about 10° may be written as 

p x q i 2 = G ( M l)P t 


(4) 


where for constant density from the shock to the stagnation point 


G(Mi) 




( 5 ) 


This function is compared with the exact solution for G^Mj) in figure 1 where it is seen 
that for Mi > 2 both relations give small changes in G for large changes in Mach num 
ber. In the following derivation the fluctuating part of G is therefore neglected. It is 
important to note that fluctuations in p are not as yet neglected except through their 
effect on the G function which is nearly constant for large changes in flow quantities. 

Equation (4) is now applied to fluctuating turbulent flow, where in general the mag- 
nitude of q^ is not the same as uj_, by assuming that the equation is correct for the 
instantaneous flow quantities. Substituting mean and fluctuating quantities for these 
instantaneous values in equation (4) and taking the time mean then gives 




p T u T 

pu 


( 6 ) 


where it has been assumed that v = w = 0 so that q 2 = + 2uu T + q’ 2 . Note that the 

subscript 1 has been dropped and that third- and higher order products of fluctuating 
quantities are neglected in equation ( 6 ) and in the following equations. If a suitable rela- 
tion between u r2 and q r2 could be provided, equation ( 6 ) could be solved for j/u 12 by 
using a correlation-coefficient form for p f u T . However, the dependence on q r2 , as well 
as on the mean flow quantity Gp t /pu 2 , can be eliminated by utilizing the additional experi 
mental data for pj 2 . That is, by first squaring both sides of equation (4), substituting 
mean and fluctuating quantities, and finally taking the time mean, the following equation 
is obtained: 


°V t 


p^u^ 


V 



1 +4 



+ 2 




( 7 ) 


I G ^tf 

The quantity I — 5 - is then eliminated by the use of equation ( 6 ), and the second-order 

\pu z / 

correlation coefficient 



is introduced. The resulting expression for the intensity of the fluctuating velocity may 
be written as 


U = 


_ -Ar + VA^r^ + - r^ 

2 


(9) 


where only second-order products of fluctuation quantities are retained. It may be noted 
that U is a function only of P and r (the ratios of mean-square fluctuations to mean 
values of pitot pressure and density, respectively) and of the correlation coefficient A. 
The effect on U of the likely range of values for the coefficient A is evaluated in 
appendix B. As mentioned previously, the magnitude of the fluctuating vector q' does 
not appear in this final equation. Also, a relation identical to that of equation (9) is 
obtained if qj is replaced by u-^ in equation (4) and if all terms of third and higher 
order are consistently neglected (see ref. 14). It can therefore be concluded that to the 
order of approximation used in equation (9), the fluctuations in mass flow caused by the 
cross velocity components v’ and w' have no effect on P in spite of the fact that, 
in general, the magnitude of the instantaneous vector q is not equal to u. 

If third-order terms are retained in the derivation, certain third-order correla- 
tions involving q' , u', and p' are present in the final equation for U. The effect of 
these third-order correlations on U is evaluated in appendix C. It is shown in appen- 
dix C that over the maximum possible range for these third-order correlations, a maxi- 
mum uncertainty of about ±20 percent in predicted values of U is possible. 


A second and more simplified method to derive an approximate relation between 
the fluctuating velocity vector, density, and pitot pressure is to take the derivative of the 
log of equation (4) (with q2 = u^ + v^ + w^) which gives 


d£ + 2 H_dH + 2 + 2 

P q2 q2 


w dw _ dG ^t 

q2 ‘ G P t 


( 10 ) 


Then, as in reference 3, for example, the differential quantities are considered the same 
as the corresponding fluctuating quantities, and all other quantities are considered as 
mean values. After squaring both sides and taking the time mean, equation (10) may be 
written as 


r p u 




(ii) 


After the introduction of the coefficient A, the solution of equation (11) for U is identi- 
cal to equation (9). 


8 



The only limitations to the generality of equation (9) are that fluctuation terms of 

G’ 

third and higher order are neglected and that -=• « 1.0. From figure 1, this latter 

G 

assumption restricts the use of equation (9) to conditions where > 2.0; that is, the 
equation can be used over all the hypersonic turbulent boundary layer except near the 
wall < 2). The data of reference 6 indicate that the local Mach number is about 2.5 
at a value of y/5 of approximately 0.05. Therefore application of equation (9) to the 
data of reference 6 is expected to be valid for ^ § 0.05. 

RESULTS AND DISCUSSION 

Measurements of Density and Pitot-Pressure Fluctuations at Mach 9 

Figures 2 and 3(a) show the variations in the ratios of root-mean-square values of 
the density and pitot -pressure fluctuations (also shown in ref. 6) to the respective local 
mean values across the boundary layer. Figure 3(b) shows the ratio of the root-mean- 
square value of wall -pressure fluctuations to the mean dynamic pressure at the boundary- 
layer edge. The data in figures 2 and 3(a) and the data from reference 6 in figure 3(b) 
were computed at the Langley Research Center from the original oscilloscope records of 
the density and strip-chart records of pitot and wall pressures. 

Under the assumption that the sensors used in reference 6 have a linear output 
over the range of a fluctuating signal, the ratios of the root-mean-square values of den- 
sity and pitot -pres sure fluctuations to the respective mean values are essentially inde- 
pendent of the absolute levels of density or pitot pressure. Consequently, the accuracy 
of the individual data points in figures 2 and 3(a) is believed to be better than the accu- 
racy of the mean data as given in reference 6. The fluctuations of both density and 
pitot pressure increase rapidly near the wall. The density profiles near the wall indi- 
cate a definite peak which, by analogy with incompressible data, is presumably in the 
vicinity of the edge of the sublayer. The overall level of the pitot-pressure fluctuation 
is considerably higher than that of the density fluctuation because, as shown by equa- 
tion (11), the fluctuating pressure depends essentially on the sum of the fluctuating den- 
sity and velocity and on the correlation p'u', which is always positive as indicated in 
appendix B. 

It should also be noted from figures 2 and 3(a) that a definite minimum in the dis- 
tribution of the fluctuation intensities of both density and pitot pressure generally occurs 
at values of y/6 from about 0.2 to 0.5. Shown in figure 4(a) are the data points and 
faired distribution of the ratio of mean density to the density at the edge of the boundary 
layer as previously presented in reference 6. Presented in figure 4(b) are the same data 
but with a possible alternate fairing which is used herein. If it is assumed that the den- 
sity fluctuations can be modeled by a mixing-length expression as previously discussed, 
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then the "flattening out" in the faired mean density profiles shown in figure 4(b) would 
imply a minimum in the density fluctuation profiles at about the same locations as shown 
in figure 2. It is of further interest to note that similar "flat" regions in static tempera- 
ture profiles have previously been reported (ref. 15) and have been observed in data for 
a turbulent boundary layer on the wall of a Mach 8 tunnel at the Langley Research Center. 


The intensity of fluctuation in the wall static pressure /p^y obtained from 
original records is plotted in figure 3(a) (see solid symbols). If the limiting value of 


obtained from the 


yPf / p t y = 0 this quantity, then the data of figure 3(a) indicate that the peaks in 
the pitot -pressure fluctuations are very sharp and close to the wall. These wall-static- 
pressure data were obtained with the same type of transducer used for the pitot -pressure 
measurements; however, attenuation of the fluctuating wall pressure may have occurred 
to a greater extent than that of the pitot pressure. This greater attenuation might be 
expected because of the lower absolute level in static pressure and differences in volume 
and lengths of tubing connecting the transducers to the pressure orifices. 

The fluctuations in wall static pressure are compared with previous data in fig- 
ure 3(b) where ]jp ^ P e u e 2 is plotted as a function of M e . The data in this figure are 


is plotted as a function of M e . The data in this figure are 


from a number of different sources and are shown here as presented in figure 4 of refer- 
ence 16 except that the data from reference 6 and the line attributed to Laufer (ref. 17) 
have been added. This line was obtained directly from lines faired through the mean of 

the data of figures 2 and 3 in reference 17; j^p^/4- p g u e ^ was obtained from figure 2 and 

/Pw 2 // p e 2 f rom figure 3 of that reference. The present data (from ref. 6) are lower than 
all other data; however, the indicated trend with Mach number is consistent with the avail- 
able data. 

As a further Chech on the level of data for JS the correlation elation 

Y W 

p s 

/ ■ w 9 = 0.006 ^ (12) 

in u Z *e 

2 e 


may be derived from the equation \J « 0.006 i pu e ^ for incompressible flow (from 
ref. 18, given in notation of present paper) by replacing the constant density p with the 
mean density at the edge of the sublayer. Equation (12) has been evaluated for adiabatic 
wall temperatures and for the cold wall temperatures of reference 6, and the results are 
shown in figure 3(b). It is seen that most of the previous data which were for adiabatic 
wall conditions are below the adiabatic -wall correlation-equation line and that the data 
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of reference 6 are considerably below the corresponding cold-wall correlation-equation 
values. Values of P s JP e used in these cold-wall estimates were taken as 0.35 to 0.6 
from figure 4(b). 


Evaluation of Mixing-Length Approach 


A comparison of the root-mean-square density fluctuations computed by equa- 


tion (1) |where the values of 


(p/pe) were taken f rQ m fig. 4(b)] with the experimental 
d(y/6) / 

data of reference 6 is shown in figure 5. The bands include all data from the four test 
runs indicated in figure 2. In general, the mixing-length approach predicts the intensity 
of the density fluctuations to within about a factor of 2. 


The root-mean-square temperature fluctuations of Kovasznay (ref. 2) and of 
Kistler (ref. 3) have also been compared with calculations obtained by a mixing-length 
approach. Calculations of the fluctuating temperature data, shown in figure 6(a) for ref- 
erence 2 and figure 6(b) for reference 3, were made by using the slopes of the mean tem- 
perature data and equation (2). The total temperature through the boundary layer was 


taken to be constant for reference 2 in forming the fluctuation temperature ratio i/t^/t, 
(see fig. 6(a)). The mean total temperature and the mean velocity distributions across 
the boundary layer were used to reduce the data of reference 3 shown in figure 6(b) for 
1.72 g M e g 4.67. For both sets of data the mixing-length approach generally predicted 
the intensity of the temperature fluctuations to within somewhat better than a factor of 2. 


A comparison of the root-mean-square velocity fluctuations computed by equa- 
tion (3) with those computed by equation (9) for the data of reference 6 are shown in fig- 
ure 7. Both hatched bands include all computed points from the four test runs. Velocity 
derivatives were obtained from fairings of mean-velocity profiles shown in figure 8. 
These velocity profiles were computed from the faired density profiles of figure 4(b) and 
the original data for in reference 6. Results from equation (9) were obtained from 
the fluctuation intensity distributions of figures 2 and 3 with A = 1.0 and with u/u e 
from figure 8. These results from equation (9) are discussed in more detail in the next 
section. The mixing-length approach (eq. (3)) predicts velocity fluctuation intensities 
that are roughly a factor of 3 larger than the results calculated from equation (9). 

The root-mean-square velocity fluctuations of Kov&sznay (ref. 2) and Kistler 
(ref. 3) have also been computed from the mixing-length equaiion (eq. (3)). The results 
of these computations are compared with measured values in figure 9 where the mean 
velocity data through the boundary layer were taken from the respective references. 

The bands for the reference 3 results (see fig. 9) are representative of three tests 
(1.72 g M e = 4.67). For both sets of data (refs. 2 and 3), the mixing-length calculations 
again are correct to within a factor of 2. 
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Correlations of Root-Mean-Square Velocity Fluctuation 


Distributions across the boundary layer.- Figure 10 shows the longitudinal velocity 
fluctuation intensities determined from equation (9) by using the distributions of fluctua- 
tion quantities in figures 2 and 3 and A = 1.0. In appendix B it is shown that the coeffi- 
cient A should be positive and values of A would be expected to ra nge from about 0.5 

i . r\ i i _ j ' * 1 • t~» • i • j _ „ : „ -T r „ ,, i l/T»2 /— — r , i n j 


to 1.0; calculations in appendix B indicate a m axim um effect on Vu '/u of ±10 percent 

for this range of A values. The ratios of l/u ,2 /u e (fig. 10) were obtained by utilizing 
the u/u e values shown in figure 8. Also shown for comparison in figure 10 are data 
obtained by hot-wire techniques (refs. 3, 19, and 20). The line shown in figure 10 for 
reference 3 represents a mean of data at three different Mach numbers; only one line is 
shown because the data were correlated to within a relatively narrow band by the use of 
u e rather than u T as the normalizing parameter. To illustrate this improved correla- 
tion, the data of reference 3 normalized with respect to u T are plotted in figure 11(a) 
as originally presented in reference 3, and the same data normalized with respect to u e 
are plotted for comparison in figure 11(b). The data are correlated to within about 
30 percent over the entire boundary layer by the use of u e . This improved correlation 
indicates that the root-mean-square velocity fluctuations normalized by the external 
velocity may be independent of Mach number and should probably be used instead of the 
local energy of longitudinal fluctuations normalized by the wall shear stress as proposed 
by Morkovin in reference 1. 

In figure 10 it can be seen that near the wall the velocity fluctuation levels for the 
present data from reference 6 (from eq. (9)) show a large increase to a peak value that 
is presumably in the vicinity of the sublayer edge and that appears to be somewhat farther 
from the wall than the peak values for the low-speed data of references 19 and 20. In the 
outer portion of the boundary layer > 0.6^ the present results are in reasonable agree- 
ment with previous data. The theory of reference 13 predicts an inverse relationship 

between {u^ju e and Rg; thus it might be expected that the results from reference 6 
should be at a higher level (for the same y/ 6) than the data of references 3 and 20 since 
values of Rg for the data of reference 6 are smaller. However, the present results in 
the range 0.2 = ^ = 0-6 are lower than previous data, either because of a Mach number 
or wall-temperature effect or because of some anomaly in the experimental conditions of 
reference 6. The trends shown in figure 10 are again largely dependent upon the trends 
of the fluctuating data in figures 2 and 3(a), which were used to compute the fluctuating 
velocity. If the trends shown for these fluctuating data are realistic, then possible expla- 
nations for the smaller values of )/u' 2 /u e noted may be as follows: (1) The smaller 

values may result from finite flow disturbances which could affect the mean profiles 
upstream of the test section; or (2) the large favorable pressure gradients upstream 
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of the test section may have reduced the value of yu'^/u e . That is, the midregion of 
the boundary-layer profiles from reference 6 may not be representative of flat-plate 
turbulent-boundary-layer flows. If the data in this region are considered atypical, the 
general agreement elsewhere between the present Mach 9 results and the previous lower 
Mach number data (fig. 10) indicates that Morkovin's hypothesis in reference 1 (that the 
basic turbulent mechanisms are independent of Mach number) can provisionally be 
extended to the lower range of hypersonic Mach numbers and cold wall conditions. 

Variation with R eyno ld s number .- In comparing the available longitudinal velocity 
fluctuation data at different Mach numbers in the present report, an attempt was made to 
use data taken at similar values of Rg so that the effect of Mach number could be 
determined independent of any Reynolds number effects. It would also be of interest to 
examine the available velocity fluctuation data, most of which is at subsonic speeds, for 
any effect of Rg, both for use in a comparison with the prediction method of reference 13 
and for possible later extrapolation to higher speeds. 

The results obtained from a survey of the available data (from refs. 2, 3, 6, and 19 
to 38) are shown in figure 12. Values of the Mach number are inserted next to the data 
points when the Mach number is appreciably different from zero. Shown for comparison 
with the data are predictions of reference 13 for incompressible flow on a flat plate. 

Since the computed quantity in reference 13 is total turbulent kinetic energy, some rela- 
tion between )/u'^ and the other turbulent velocity components must be postulated before 

a prediction of \Zu'2 can be made. For the near-wall region, the data of reference 19 
indicate that a reasonable relation is q'^ = 1.5u'^, which gives the das hed lin es 
(figs. 12(a) and 12(b)). The solid line (fig. 12(c)) is from the relation q’^ = 2u'^ which, 
as indicated in reference 13, is a reasonable approximation for the outer region of an 
incompressible turbulent boundary layer. 

Peak values of the velocity fluctuation data normalized with respect to u T (where 
u T values could be obtained with certainty) are plotted as functions of Rg in fig- 
ure 12(a). Most of the data are within about 20 percent of the level predicted in refer- 
ence 13 except for the present data, the data from references 34 and 38, and some data 
from reference 25 (occurring between 5 X 10^ ~ Rg ~ 2 X 10 4 ) which are believed to have 
been obtained in transitional flow. The present Mach 9 results from reference 6 
(5 X 10 4 £ Rg £ 8 X 10 4 ) fall on the lower side of the data scatter. These data points 

were obtained from the peak values of /u' ^ as computed with equation (9) for the runs 
of reference 6 for which surface shear data were available. It is seen that in spite of 
the high Mach numbers and cold wall conditions, these data are in reasonable agreement 
with previous data which are for adiabatic wall conditions and a maximum Mach number 
of about 0.7. It should be noted, however, that the viscous sublayer is evidently much 
thicker for these low Reynolds number, Mach 9 data than for the low-speed data. That 
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is, if the thickness of the viscous sublayer is taken as ~ 0.1, which corresponds to 

the approximate location of the peak in l/^2/ue for the present data (see fig. 10), the 

P — u ^ r 

y s range from about 60 to 120. The cor- 


values of a law-of-the-wall coordinate 


P u_ 
p w t 

Mw 


responding range of values in this coordinate for the low-speed data of figure 12(a) is 
from 15 to 20. When the effects of variable density are included by using the Baronti- 
Libby transformation theory (ref. 39), the values of the transformed law-of-the-wall 


coordinate for the present data at 


1 ) 

6Js 


0.1 are from about 40 to 70, again a range 
which indicates an unusually thick sublayer. 

Velocity fluctuation data are usually normalized by u T rather than by u e , and the 


ratio Vu ,,a /u T has been found to be relatively invariant with Reynolds number at low 
speeds, as discussed for instance in reference 40 and shown in figure 12(a). However, 
figures 11(a) and 11(b) have indicated that the apparent effect of Mach number is reduced 
if the velocity fluctuation data are normalized by u e rather than by u r , and therefore 


i(u'2/u e may be a useful parameter for comparing low-speed data with data from com- 
pressible flows. The root-mean-square velocity fluctuation data at the peak and at 
Z = 0.5 are plotted in this form in figures 12(b) and 12(c), respectively. Correlation of 
the peak values may be slightly better when u T is used (fig. 12(a)) than when u e is 



peak 


with R5 is in 


At ^ = 0.5 (fig. 12(c)) the data are not as well correlated as 


used (fig. 12(b)). In the latter plot the predicted trend of 

agreement with the data, 
in the previous figures and the Mach 9 data are considerably lower than most of the pre- 
vious data. The trend of the theoretical prediction (ref. 13) is supported by the data but 
the predicted level is too high. There appears to be no consistent variation with Mach 
number in figure 12(c); however, the Mach 9 results are low. 

In view of the previous discussion concerning the detailed comparison of the 

^i/u e profile for M ~ 9 with lower Mach number data (fig. 10), the fact that the data 

for M ~ 9 at £ = 0.5 are low (fig. 12(c)) probably does not invalidate the tentative 
0 

extension of the hypothesis of reference 1 to hypersonic Mach numbers of about 9 for cold 
wall conditions and boundary layers with small pressure gradients. 


Correlations of Root -Me an -Square Temperature or Density Fluctuations 

The density fluctuation data from reference 6 (as presented herein) may be com- 
pared directly with temperature fluctuation data from other sources since, from the equa- 
tion of state, 

dp _ dT dp 
P " T P 
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Then if the level of pressure fluctuations is neglected (see appendix B), the root-mean- 
square ratios of density and temperature are equal 



P T 


In comparisons of /p'2 with previous data for |/t'^, a normalizing parameter 
that would minimize variations caused by Mach number and wall temperature should be 
used. For velocity fluctuations it has been shown (figs. 10 and 11) that the appropriate 
normalizing parameter is u e , which is the mean velocity at the edge of the boundary 
layer (i.e., the difference between the velocity at the wall and the edge). The same rea- 
soning applied to the temperature or density fluctuations indicates that the normalizing 
parameter should be T w - T e or p g - p w for situations where the wall temperature 
or density is the maximum or minimum value, respectively, across the boundary layer 
such as for adiabatic or hot wall conditions. For cold wall conditions where heat is 
transferred to the wall from the boundary layer, there is a peak in temperature T s and 
a corresponding minimum in density p„ very close to the wall. The normalizing param- 
eter for these conditions should then be T s - T e or p g - p g . Similar normalizing 
parameters (maximum velocity and density defects) were found to correlate velocity and 
density fluctuation data in a compressible wake flow (ref. 41). 

To determine the effect of increasing Mach number for these normalizing param- 
eters, the temperature fluctuation data of references 2,3, and 26 have been normalized 
by the difference between their respective wall and boundary-layer-edge temperatures, 
and the results are plotted as functions of y/6 in figure 13 (note that T w ~ T s for 
ref. 26 because M « o). Similarly, the fluctuation density data of reference 6 have been 
normalized by the differences between the edge and minimum density and are also shown 
in figure 13 (the ideal gas values of p g from ref. 6 have been used). The dimension- 
less fluctuating density and temperature profiles shown from the various sources are 
in agreement within a factor of about 2 over most of the boundary layer. Consequently, 
the nominal trends and magnitudes shown in figure 13 are useful for estimates of fluctua- 
tion levels in both compressible and incompressible boundary layers. 

CONCLUDING REMARKS 

Estimates of the variation in intensity of the longitudinal velocity fluctuations across 
a hypersonic turbulent boundary layer at Mach numbers from 8.2 to 8.9 have been made by 
utilizing experimental measurements of fluctuating density and pitot pressure. Analysis 
of the momentum relation used to make these estimates indicates that when all third 
and higher order correlations of fluctuating quantities are neglected, the fluctuations in 
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mass flow normal to the mean velocity vector have no effect on fluctuations in pitot 
pressure. 

Comparison of available data indicates that for the outer part of constant -pres sure 
boundary layers, the root-mean-square values of the fluctuating velocity normalized by 
the mean velocity at the edge of the boundary layer may be independent of wall tempera- 
ture and Mach number up to Mach 9. In the near -wall region where the intensity of 
velocity fluctuations reaches a maximum, the shear-stress velocity may be a more suit- 
able normalizing parameter. 

When the root-mean-square values of the density and temperature fluctuations 
across the boundary layer were normalized by the difference between the respective 
minimum or maximum mean value across the boundary layer and the mean value at the 
boundary-layer edge, available data were again found to be essentially independent of 
Mach number and wall temperature. 

On the basis of these preliminary results, it appears that by the use of mixing- 
length concepts or appropriate correlating parameters, estimates can be made of trends 
and magnitudes, within a factor of 2, of fluctuating velocity, density, and temperature for 
Mach numbers up to about 9 and for cold wall conditions. The validity of the hypothesis 
that the basic turbulent mechanisms are independent of Mach number is thereby provi- 
sionally extended to cold wall conditions and to hypersonic Mach numbers of about 9. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Langley Station, Hampton, Va., August 1, 1969. 
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APPENDIX A 


DERIVATION OF INSTANTANEOUS PITOT -PRESSURE EQUATIONS 


An approximate relation between the fluctuating velocity, density, and pitot pressure 
may be derived by assuming that the usual steady-flow equations apply to the instanta- 
neous pitot pressure. In this derivation the molecular shear is neglected and, for sim- 
plicity, the instantaneous density from the postshock point (point 2) to the stagnation point 
is assumed constant. For a conventional pitot tube at zero angle of attack as illustrated in 
the following sketch, 



the instantaneous pitot pressure measured by the tube would be (as obtained from ref. 42, 
p. 243, for example, and given in the present notation) 


p t = P lUl 2 



p l 

Pl u l 2 


1 Pl\ 

2P 2 J 


(Al) 


If the pitot tube is placed at some angle of attack such that the instantaneous veloc- 
ity vector ahead of the shock is q^ and if the shock is assumed normal to q^ as shown 
in the following sketch, 
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the measured pitot pressure will have the same dependence on the density and on the 
magnitude of the velocity vector as in equation (Al) since the entropy of the stagnation 
streamline is identical if q^ = uj. The following derivation (analogous to that of 
eq. (Al)) is included to illustrate this effect of angle of attack on pitot pressure. 

Conservation of momentum across the shock requires that 

Pi + P-fii 2 = P 2 + P 2 <1 2 2 (A2) 


The momentum equation along a postshock streamline is 


pq 



(A3) 


If the streamline shown in the preceding sketch is assumed to stagnate on the tube face at 
point t, then the measured pitot pressure is approximately (by integration of eq. (A3)) 


Pt = P2 


+ 


* 2*2 

2 


(A4) 


where again the instantaneous density from point 2 to point t is assumed constant. By 
eliminating p 2 from equations (A2) and (A4) and using the continuity equation 
Pl^l ~ *2^2’ f°ll° w i n g equation is obtained: 


P t = P^ 2 



Pi 

P^l 2 



(AS) 


Comparison of equations (Al) and (A5) shows that these equations are identical in terms 
of the local velocity vectors uj and q^, as would be required by the assumption that 
the pitot -tube shock is normal to the stagnation streamline. The measured pitot pres- 
sure would therefore be insensitive to angle of attack as long as this assumption is valid 
and the magnitude of q^ is the same as u-^. Experimental results of reference 43, for 
example, indicate that these assumptions are correct for pitot -tube angles of attack up to 
about 20° at Mach numbers up to 6, since the measured total pressure was essentially 
independent of angle of attack over this range of conditions. 

To determine the effect of pitot -tube angle of attack at higher Mach numbers, some 
data were obtained at a Mach number of 19 in the Langley hypersonic nitrogen tunnel. 
These data are shown in figure 14 where the measured pitot pressures (normalized to 
the values at a - 0) are plotted against angle of attack of the tube. (The experimental 
conditions for these data are given in the figure.) These data indicate that the effect 
of pitot -tube angle of attack on total pressure is about 7.5 percent at a = 20° for 
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M e = 19. The increased effect of a at Mach 19, as compared with the effect at lower 
Mach numbers in reference 43, may not be entirely due to the increase in Mach number, 
since the pitot -tube lip shape was different for the two sets of data. The sharp-lip con- 
figuration of reference 43 is known to minimize the angle -of-attack effect, as indicated 
in reference 44. The blunt tip used for the Mach 19 data is more nearly like the pitot- 
tube tip used in reference 6. The data of figure 14 are therefore thought to more nearly 
represent the effect of a on p^. for the conditions of reference 6. 
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APPENDIX B 


EVALUATION OF SECOND-ORDER CORRELATION COEFFICIENT A 


The purpose of this appendix is to show that definite limits on the values of the 
second-order correlation of fluctuation density and velocity (eq. (8)) can be established 
from the measurements of reference 6 f or mean and fluctuating pitot pressure and den- 
sity. It is shown that the values of u’ 2 are not overly sensitive to values of the cor- 
responding correlation coefficient as prescribed by these limits. 


Preliminary Remarks 


Morkovin (ref. 1) has shown, for boundary layers with zero heat transfer, the fol- 
lowing temperature -velocity correlation coefficient (given here in the notation of the pres- 
ent paper): 


Rtu = 



- 1.0 


This estimate has been verified to within about 50 percent by data of Morkovin at 
M e -1.75 (ref. 1) where -0.9 < Rp u < -0.7 and by data of Kistler (ref. 3) where 
-0.8 < Rp u < -0.5 up to Mach number 4.67. From the equation of state with negligible 
pressure fluctuations, the fluctuations in p and T are related by the expression 

9L= - lL 

p T 


By using this relation, it can be inferred that for adiabatic boundary layers, the values of 


the density-velocity correlation coefficient A = 


p'u’ 


would be in the same range 


( p r2 u .2) l/2 

as R-pu but of opposite sign. Although there are no direct measurements of static 


pressure 


fluctuations in a boundary layer, it is reasonable to neglect { P’Vp, which is 


'2 ?, 


p, on the basis of the low levels of VP^/p w shown in fig- 


small compared with yp 

ure 3(a) and the expectation that \j p’ 2 would decrease by an order of magnitude from the 

wall value over most of the boundary layer. This rapid decrease in ]/p^ with increasing 
y in the near-wall region is consistent with the theory of reference 45 and with the data 

of Laufer (ref. 17), who found that free-stream values of /p' 2 were one-tenth of the wall 
values. 
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The magnitude of A can be independently determined, since it can be reasoned as 
in reference 1 that for a region (even in a cold-wall boundary layer) where both 9p/ 9y 
and 8u/8y are positive, the correlation p’u' should be positive. That is, if a lump of 
gas with a smaller than average density came from a layer close to the wall, it would 
have a velocity smaller than average. Since the data of Wallace (ref. 6) satisfy these 
conditions on the mean flow gradients for most of the boundary layer as shown by fig- 
ures 4 and 8, the values of A should be within the limits 0 § A S 1.0. 


Basic Equations 

Four independent relations between the instantaneous flow quantities are required 
in order to evaluate the coefficient A. The pitot -pressure relation (eq. (4)) has been 
derived in appendix A and its limitations discussed in the text. This equation may be 
written as 

pq 2 = p(u 2 + v 2 + w 2 ) = G(M)p t (Bl) 


In place of the "Strong Reynolds Analogy" assumption used by Morkovin (ref. 1) for 
adiabatic walls, it is assumed herein that the relation 



(B2) 


represents a solution to the energy equation that accounts roughly for the effects of 
upstream pressure gradient_and of Prandtl number not equal to 1.0. Equation (B2) should 
probably be used only for > 0.5 and cold wall conditions. Values of n that might 
apply to instantaneous flow quantities are unknown, so for the present purposes it is fur- 
ther assumed that values of n corresponding to mean flow data can be used. That is, 
depending on the type of flow under consideration, n would generally be in the range 

1.0 < n < 2.5 


The other relations required are the equation of state 

p = pRT (B3) 

and the definition of local stagnation temperature 


T t = T + 


u^ 


+ v 2 + w^ 

2cp 


(B4) 


It has been shown in the text (compare eqs. (10) and (11) with eq. (9)) that second- 
order correlation equations can be obtained directly from differential forms of the basic 
relations. Before applying this procedure to equations (Bl) to (B4), it should be noted 
that multiplication of equation (B4) by pR and the use of equation (B3) give 
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p + p( u 2 + v 2 + w 2 ) = RpTt 


By taking the log derivative of this equation and rearranging terms, the equation 

^t = d£-d£ + ( y _ 1)M 2 ^ 

TP P u 


(B5) 


is obtained when v = w = 0. Similarly from equation (B2) (for T w assumed constant) 




(B6) 


and from equation (Bl) 


dp 


i= d£ + 2 *£ 

n U 


(B7) 


Correlation Equations 


In the application of equations (B5), (B6), and (B7) to turbulent flow, the restriction 
is that the maximum deviations from the mean are small, since the derivatives are 
replaced by the corresponding turbulent fluctuation quantities. These equations can, 
nevertheless, be manipulated in various ways depending on what data are available and 
which terms are considered negligible. The following development is for the present 
situation where data for p’ and p£ are available and p’/P or P' 2 /P 2 are generally 
negligible compared with other terms. 

By first neglecting the pressure term dp/p in equation (B5), squaring the resulting 
equation, and taking the time mean, the following expression is obtained: 



+ a 



(B8) 


2 

In this equation a = (y - 1)M . Equating this result to the square of equation (B6) then 
gives 


£l-2a£^ + (a 2 -b 2 )^i=0 

p 2 P u u 2 


(B9) 


where b = n^l + - - - 

tion (B7), another relation between 


. By squaring and taking the time mean of equa- 
p' and u' is obtained as 
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u' 2 _lfPt P_\ W 

S 2 "% 2 p 2 ~P 5 


(BIO) 


Equations (B9) and (BIO) are then solved for the two unknowns p'u'/pu and u' 2 /u^ to 
give the following expressions: 


The coefficient 


p’ 2 . 1/ Q h \2j Pf p^ I 

P \Pt P 


pu a 2 

' - b 2 h 

u' 2 _ 1 


u 2 a 2 - b 2 + 2a 

V 

I 


It__pll.pl 


-2 

P 


(Bll) 


(B12) 


A = 


p'u' _ p'u’ 


pu 


u 


can then be evaluated and, as an indication of the reliability and consistency of the data, 
the values of A should be within the limits 0 g A s 1.0 which were mentioned previ- 


ously. Examination of equation (B12) shows that for 


x w 

Tt 


« 1.0, corresponding to very 


cold wall conditions, negative values of u’^/u 2 are predicted for n g 2.0. For the con- 
ditions of reference 6 where the minimum T w /Tj. was about 0.09, it therefore seems 
reasonable to restrict n to the range 1.0 < n < 2.1. 


Application to Data 

Values of the correlation coefficient A for the four runs of reference 6 as used 
herein have been computed from equations (Bll) and (B12) for n = 1.6 and n = 1.8. All 
these data gave values in the range 0.4 < A < 4.5 for ^ > 0.1. Since values of A > 1.0 
are not physically possible, it can be assumed that data scatter is responsible for these 
large values. The equations thus predict that for these values of n, 0.4 < A < 1.0. It is 
therefore justifiable to assume A is constant, whereupon equation (B10) becomes quad- 
ratic in P/u and equation (9) of the text is obtained. Values of A in this range have 


been used to compute yu’ 2 / u from equation (9) for one of the runs of reference 6, and 


23 



APPENDIX B 


the results are shown in figure 15. It can be seen 
dieted values of l/u'^/u over this range of values 


that the maximum effect on the pre- 
for A is about ±0.01. 
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POSSIBLE EFFECTS OF THIRD-ORDER CORRELATIONS 


After a reasonable range of values is established for the second-order correlation 
coefficient A, then possible effects of third-order correlations on the values of the pre- 
dicted root-mean-square velocity fluctuations may be examined. 


If the same procedure used in the text to derive equation (9) is followed, except that 
all terms consisting of up to third-order products of fluctuation quantities are retained, 
the following equation is obtained: 


u' 


-2 

u 


u’q' 2 [ 2 p'u’ 2 


u 3 


pu 2 



p ,2 u' 
P U 




p'u* 

pu 


(Cl) 


Comparison of this result with equation (BIO) shows that four third-order terms that 
involve u', p' , and the fluctuating velocity vector q' = ( u '2 + v '2 + w ' 2)*/ 2 h ave been 


added. In order to write an explicit expression for 
coefficients are defined as 



the third-order correlation 



R 

pu 


2 


R pq 2 


p’u’ 2 
{p^ u’ 2 

p'q' 2 _ p’q' 2 

\fp' 2 q'2 c \/ p '2 u’ 2 


R 2 

p^u 



J 


(C2) 


where 



Introduction of these coefficients into equation (Cl) gives the following 


cubic equation for U: 


C V uS - l 1 + 2 V r + i R pq2 r ) u2 + ( Ar + R p2u r2 ) U - ?( p2 - r2 ) = ° < C3 > 
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Since the sign and exact magnitude of the third-order coefficients are unknown, the 
best that can be done at present is to obtain solutions to equation (C3) for a range of con- 
stant values assigned to the third-order correlation coefficients. These solutions have 
been carried out with A = 1.0 and c = 2 (as used in ref. 13). The range of values 
used for the four third-order correlation coefficients was varied independently from -1.0 


to 1.0 corresponding to perfect anticorrelation and perfect correlation. The results are 


shown for a typical run in figure 15. When 


the effect of these coefficients is 


maximized at about ±20 percent. When 



~ 0.02, the effect is less, as would be 


expected, and amounts to about ±10 percent. The predicted values of Vu’2/u were the 
largest when all coefficients were -1.0 and the smallest when all coefficients were +1.0. 
Since these particular combinations of limiting values are not likely, it can be concluded 
that unless the intensity of velocity fluctuations considerably exceeds 0.1, the third-order 
correlations are indeed negligible. 
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Figure 2 - Distribution of measured density fluctuations. 
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(a) Measured pitot-pressure fluctuations. 
Figure 3.- Pressure-fluctuation distributions. 
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(b) Comparison of wall-pressure fluctuations in turbulent boundary layer. 
Figure 3.- Concluded. 
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(b) Present fairing. 

Figure 4.- Mean-density profiles. 


35 



\ 0 \ 







y/S 

(b) Comparison with data from reference 3. 
Figure 6.- Concluded. 
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Figure 7.- Comparison of present fluctuating velocities computed from equations (3) and (9). 
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Figure 8.- Mean-velocity profiles. 
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Figure 9.- Comparison of measured fluctuating velocity from previous experiments (refs. 2 and 3) with mixing-length calculations (eg. (3)) 
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Figure 10.- Comparison of present velocity fluctuation intensity with previous results. 
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(a) Normalized by shear-stress velocity. 

Figure 11.- Fluctuating velocity through turbulent boundary layer (taken from ref. 3). 
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(b) Normalized by u e . 
Figure 11.- Concluded. 
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(c) Velocity fluctuation (normalized by Ue) at | = 0.5. 
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Figure 12.- Concluded. 
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Figure 14.- Effect of angle of attack on pitot pressure. 
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